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ABSTRACT 

A finite element Galerkln formulation has been developed to 
study electromagnetic propagation In complex two-dimensional 
absorbing ducts. The reflection and transmission at entrance and 
exit boundaries are determined by coupling the finite element 
solutions at the entrance and exit to the eigenfunctions of an 
Infinite uniform perfect conducting duct. Example solutions are 
presented for electromagnetic propagation with absorbing duct 
walls and propagation through dielectric-metallic matrix 
materials . 


INTR0DUC1 ION 

For Inlets and exhaust ducts of jet engines, the structural 
design and materials are usually determined by aerodynamic, ther 
modynamlc, and structural requirements. In recent military 
applications; however, electromagnetic constraints have begun to 
play a major role In determining the structural composition and 
geometrical configurations of Inlet nacelles and exhaust ducts. 
Clearly, the aerodynamic and electromagnetic disciplines must 
Interact In an effective manner to produce an efficient propul 
slon system. 

The NASA Lewis grant program In electromagnetics (Refs. T 
and 2) has been concerned with a basic understanding of electro 
magnetic propagation In Internal propulsion systems. Following 
the same theme, the present paper presents a finite element 
analysis potentially capable of evaluating the electromagnetic 
absorption and scattering properties of jet engines composed of 
complex structural geometries and composite materials. 

ME1H0D OF ANALYSIS 

The propagation of electromagnetic waves are studied In the 
two-dimensional (y, z) finite element network shown In Fig. 1. 
This configuration Is Ideal for modeling the grazing absorption 
by wall materials or direct penetration through materials Inside 
the duct. The calculation domain Is composed of uniform 1 nf 1 
nltely long entrance and exit regions with perfectly conducting 
walls and the central, variable property, nonuniform region. 


In the absorbing variable property region of Fig. 1, the 
finite element analysis is employed to determine the field var 
iable. For harmonic wave propagation in variable property mate 
rials, Maxwell's equations can be combined to yield a wave 
equation of the form 



where H represents the magnetic field, c the complex per- 
mittivity, v the permeability and w the wave frequency. For 
TM (transverse magnetic) wave propagation in the two dimensional 
geometry of Fig. 1, Eq. (1) reduces to 

V • ~ WH x f w 2 vH x = 0 (2) 

The boundary conditions associated with Eq. (2) require 
setting the normal gradient of H x to zero at the perfectly 
conducting outer walls and matching the finite element solutions 
at the entrance and exit planes of the nonuniformity to the well- 
known analytical eigenfunction expansions in the uniform inlet and 
outlet ducts. This permits a multimodal representation accounting 
for reflection and transmission at the ducts inlet and outlet. 

Equation (2) and the appropriate boundary conditions have 
been solved by Finite Element lheory using the Method of Weighted 
Residuals. Details of the solution procedure can be found in 
Ref. 3. 


RESULTS AND COMPARISONS 

For theory and code validation, the finite element solution 
is first applied to a case where an exact solution exits. Next, 
the finite element solution will be applied to problems with 
grazing absorption and normal Incidence on composite materials. 

The results are expressed in terms of nondimensional quantities. 
The axial coordinates have been scaled by the duct height, the 
permittivity and permeability by the free space values and the 
frequency by the duct height and the speed of light. Also, j 
represents the square root of -1. 

Example 1 Reflection and Transmission with Normal Incidence 

The first case considers a step change in dimensionless per 
mittlvity from 1 to 4 at the dimensionless axial position of z 
equals 0.2b inside the finite element grid. As shown in Fig. 2, 
for an Incident plane wave of dimensionless frequency 2ir, the 
finite element and exact analytical theories are in good agreement 
for the absolute magnitude of the magnetic Intensity H x . Ref 
erence 3 contains additional validation examples. 

Example 2 Grazing Wall Absorption 

Figure 3 shows a configuration with wall absorbers placed at 
the bottom and top of the duct for materials with a dimensionless 


2 



permeability of 4.1 and permittivity of 1.-J2.83. Again for an 
incident plane wave of dimensionless frequency 2ir, the absolute 
magnitude of the magnetic field is displayed in the duct. A five 
mode expansion in the entrance and exit ducts was sufficient to 
obtain convergence. As seen in Fig. 3, the magnetic field has 
decreased to near zero at the walls but remains relatively high 
in the center of the duct. 

Example 3 Microscopic Absorption and Reflection 

lhe previous two examples employed the macroscopic properties 
of y and c. However, many advanced coating material and 
structures employ small ferrite particles for absorption or slen 
der metallic rods for structural strength. In these cases, aver 
age values of y and e are difficult to estimate. However, 
with the large storage capacity of the modern computer, small test 
samples of the detailed material structure can now be modeled on 
the computer. Figure 4(a) displays a sample with two small iron 
squares 0.0005 in. along a side embedded in an epoxy of dimen 
sionless e = 2.0 - jO.OOl and y = 1.0. For an incident plane 
wave of 9 GHz, the magnetic field lines are shown in Fig. 4(b). 

In this case, comparison of the input plane wave to the reflected 
modes at the entrance indicate that 53.3 percent of the incoming 
energy was reflected back down the duct. Without the iron rods, 
the reflected energy due to the base material was negligible. 

This example illustrates how hidden internal metallic structural 
material can significantly effect reflection coefficients. 

CONCLUDING REMARKS 

The finite element method has been developed to handle the 
problem of electromagnetic propagation in a duct with varying wall 
properties and geometries. The numerical formulation is rela 
tively simple to use and appears to give very accurate results. 

For higher frequencies many wave lengths will exist in the duct 
requiring a large number of elements. With the expanded memory 
capability of the modern computer, the restriction does not 
represent a severe limitation of the method. 
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Figure 1.- Two dimensional wave guide finite element model. 



Figure 2. - A comparison of the magnitude of the axial 

MAGNETIC INTENSITY VARIATION ALONG THE LOWER WALL IN 
A UNIFORM DUCT WITH PERFECTLY CONDUCTING WALLS AND A 
CHANGE IN PROPERTIES OF THE MEDIA AT Z OF 0.25 (& 1 = 

1.0 AND £ 2 = 4.0) AS OBTAINED BY USING AN EXACT SOLU- 
TION AND A FINITE SOLUTION FOR A PLANE WAVE (MODE-ONE) 
INCIDENT AT Z = 0 WITH U = 27T. 



Figure 3.- Effect of upper and lower wall absorbers 

ON THE MAGNITUDE OF THE MAGNETIC FIELD FOR A FIVE- 
MODE MODAL EXPANSION IN THE ENTRANCE AND EXIT DUCTS 
(U = 27T). 
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(a) Iron rod location. 



(b) Magnetic field contours. 


Figure 4.- Magnetic field contours for a plane 9 GHz wave impinging 

ON TWO IRON RODS IMBEDDED IN A EPOXY FIELD. SAMPLE SIZE IS 
0,010 IN. BY 0.010 IN. 
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